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Abstract 

The possibility of materials that are governed by a fixed point related to small world networks is 



discussed. In particular, large-scale Monte Carlo simulations are performed on Ising ferromagnetic 
models on two different small- world networks generated from a one-dimensional spin chain. One 



has the small-world bond strengths independent of the length, and exhibits a finite-temperature 
phase transition. The other has small-world bonds built from atoms, and although there is no 
finite-temperature phase transition the system shows a slow power-law change of the effective 
critical temperature of a finite system as a function of the system size. An outline of a possible 
synthesis route for quasi small-world nanomaterials is presented. 
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I. INTRODUCTION 



Many materials have novel physical properties that are attributed to the quasi-low- 
dimensional nature of their structure, as can be seen in recent research in a number of 
areas. Examples of quasi-one-dimensional materials include conductors and superconduc- 
tors Two-dimensional thin films lead to novel effects, including the quantum Hall effect 
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and low-dimensional magnetism The change from one effective dimension to another 
may lead to interesting physical effects. For example, it may be responsible for the onset 
of high-temperature superconductivity j^]. This type of dimensional crossover can lead to 
effective non-integer dimensional fixed points |5] . However, quasi- low-dimensional materials 
have constraints due to their low-dimensional behavior. These include the absence of a phase 
transition both in one-dimensional systems and in some two-dimensional systems where two 
is below the lower critical dimension. 

Recently there has been a great deal of interest in networks that are not regular lattices, 
such as small- world networks p. The study of these networks has been motivated mainly 
by social organizations (such as six degrees of separation) and connectivities of computers 
(such as the scale-free world- wide- web network) [6]. Such connectivities have also been 
used, for example, in non-trivial parallelization of short-ranged discrete event simulations 
11]. Simulations of models such as Ising spin models on these networks have 
been studied, but no attempt has been made to ask whether these theoretical models could 
actually be designed and built via various synthesis routes. This question is addressed in 
this manuscript. There is a difference between the models studied to date and the question 
of whether or not materials can be made that are governed by small-world fixed-points. 
The difference is that materials must be constructed from atoms and must be embedded in 
three-dimensional space. 

II. MODEL AND METHODS 

The models studied here are Ising models with Nq spins on a linear chain, with a nearest- 
neighbor ferromagnetic interaction constant J\. 

In the first model, model 1, if there is a small- world connection between spins i and j 
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and i>j+la ferromagnetic interaction of strength J2(i,j) is added. The Hamiltonian is 

H = -Ji^ajCTi+x - y^ j J 2 {i,j)(Ti(Tj (1) 

i SW 

where periodic boundary conditions are used (spin N + 1 equals spin 1) and the Ising spins 
a = ±1. Terms in the second sum are only present if there exists a small world connection 
between spins % and j. We construct a (random) small- world network algorithmically by: 
i) start with spin 1, and randomly connect it to any of the other N Q — 1 spins; ii) if spin 
2 is not connected to spin 1, then randomly connect it to one of the Nq — 2 spins that 
are not already connected; iii) continue for all Nq spins. Note that here Nq must be even. 
This algorithm gives each spin i three connections, two of strength J\ to spins i — 1 and 
i + 1 and one of strength J2 to the small- world connection between spin i and j. This type 
of small-world connection has been used in the past to obtain perfectly scalable parallel 
discrete-event simulations [7]. A study of the Ising model on a similar small world network 
with J2{i,j = J r ij a (i- e - having a power-law dependence on distance) has recently been 
performed 12 1. There they conclude that any non-zero value for a destroys the finite- 
temperature phase transition in the thermodynamic limit. This result differs from the case 
with J2(i,j) = Jo (independent of length), where a finite-temperature phase transition has 
been observed |l3j, llJ, LUj lly, . Our first model has all small- world bonds of strength J2, 
independent of length. 

In our second model, model 2, we construct the small- world connections from Ising spins. 
In particular, for each small-world bond constructed as in our first model between spins 
i and j, we add + 1 Ising spins and couple each of them together with interactions of 
strength J2. These additional Ising chains are coupled to the original Ising spins i and j with 
interaction strength J2. Thus, the small-world connections have been constructed using the 
same lattice constant as the original lattice (but different nearest-neighbor interactions). 

We have performed a standard importance-sampling Monte Carlo simulation 18]. We 
measured the magnetization m = ^(Tj/iVtot, order parameter, susceptibility x, internal 
energy, specific heat, and the Binder 4 th order cumulant for the order parameter U = 

(m 4 ) 

1 — ^ 2^ ■ The code was run using trivial parallelization on a cluster, using up to 128 
processing elements. The parallel random number generator SPRNG 1.0 was used. The 
spins to be updated were chosen randomly, a Glauber spin-flip probability was used. For 
each temperature T, 10 5 Monte Carlo steps per spin (MCSS) were used for thermalization 




FIG. 1: The total number of Ising spins iVtot as a function of the number of spins Nq in the one- 
dimensional chain is shown. The solid line is a fit to the data for Nq > 2 10 . Also shown is the 
standard deviation of iVt Q t from the mean (squares), and a fit to those data (dashed line). 

and averages were taken over 10 6 MCSS (with measurements taken every 10 MCSS). The 
energy units were chosen such that Jx/ks = 1, where k-g is Boltzmann's constant. In this 
paper the choice J 2 = 4 Ji was made. 



III. DATA AND ANALYSIS 

In our second model, one question is the average number of total Ising spins N tot starting 
with a spin chain of length Nq. For a chain of length L and lattice constant a, one has 
L = Noa. Here we take the lattice spacing a — 1, so Nq — L. It is shown in Fig. 1 that 
iV to t oc Nq with exponent x ~ 2. A fit to data (averaged over 10 3 different small- world bond 
connections) for 2 10 < iVo < 2 20 gives iVtot = 0.127A^q- 999 . A fit with the same range of A^o 
for the standard deviation of the mean gives std(A" tot ) = ClOSA^- 497 . 
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FIG. 2: The internal energy per Ising spin for (a) model 1, and (b) model 2. These are parame- 
terized by the chain length, L = Nqcl, where a is the lattice spacing and Nq is the number of spins 
in the one-dimensional backbone. 

The relationship N tot oc Nq can be understood easily in an approximate fashion, by 
ignoring correlations between the small-world interconnections. Although the number of 
spins scales as the square of the linear distance Nq, the fixed point governing the system is 
not expected to be the two-dimensional fixed point since the bond connectivity is different 
from that of a two-dimensional lattice. The average length of a connection is N a/A. This 
is because the lengths of the connections are chosen uniformly up to the maximum possible 
distance of a connection, which is N a/2 due to the periodic boundary conditions. There 
are N /2 such connections (the factor of 2 is to take into account double counting). Ignoring 
factors in the number of spins that are proportional to N Q thus gives N tot pa ^Nq ■ As seen in 
the fit in Fig. 1, both the fitted exponent and prefactor are in excellent agreement with this 
argument for N > 2 7 . An extension of this argument shows that if only a fraction p C onnect of 
spins were connected, then the total number of spins should scale like N tot pa Pc ™ nect Nq for 
large Nq. Simulations were performed with values of A^ = 8 and for N up to Nq = 16384 
for model 1 and Nq = 256 for model 2. 



Figure El shows the average internal energy per Ising spin, in units of J±, as a function 
of temperature. The internal energy is the expectation value of the Hamiltonian, E = ((H). 
Model 1, Fig.|2^a), has the low-energy value for E/N independent of N, since in the ground 

5 



0.7 




0.7 



0.6 




1 0.4 

o 



0.3 - 



CD 



C 0.2 

in 



0.1 - 



o 







FIG. 3: The Binder fourth-order cumulant for (a) model 1, and (b) model 2. 



state (all spins the same direction) the energy is E/N = (2Ji + J%)/2 where the division 
by 2 takes into account double counting. For model 1, J\ — 1 and J2 = 4, so E/N = —3 
at low temperatures. Note that model 1 has a singularity developing with large N near 
T = 2.5J\. Model 2, Fig. EJb) , looks different. Note that for model 2, this energy is divided 
by the total number of Ising spins, N tot . The ground-state energy per spin has the form 
e = No(2j 1 +j2)+2MN tot -N ) which at 1qw temperatures approaches -4 since N tot > iVo for 

large N . 

Figure El shows the Binder 4 t/l -order cumulant for the magnetization as a function of the 
temperature for both models. The temperature where two of these curves for different N 
cross is an estimate for the critical temperature. An alternative estimate can be obtained 
by choosing a value such as 0.5, and letting the estimate of the critical temperature be 
where the cumulant crosses the value of 0.5. (Ideally, one would like to use the value of 
the infinite-lattice cumulant, but this is often unknown unless the universality class of the 
model is known.) Model 1, Fig. 0(a), shows a very nice indication of a critical point near 
T « 2.7Ji. Model 2, Fig. 0(b), shows no such estimate of the critical point. Only by 
choosing the temperature where the cumulant crosses 0.5 can an estimate of the finite-size 
critical point (if it exists) be obtained. Note that for model 2 there are a number of different 
quenched small-world bond configurations for N = 8, 16, and 256. 
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Figure |3] shows the susceptibility x P er Ising spin times the temperature T. For model 1, 
Fig. EKa)> shows a sharpening maximum at T « 2.7. Figure for model 1, shows the 
maximum value for \T as a function of A^ . There are five different quenched small-world 
connections shown for each value of N . The slope of this log-log plot gives an estimate 
for the exponent ratio 7/V, with 7 the susceptibility exponent and v the correlation length 
exponent. The value obtained for the slope is 0.73 using all lattice sizes, and it is 0.65 using 
only the sizes 256 and 512. Although simulations for larger N may give lower slopes yet, 
we estimate that j/u — 0.6 ± 0.1. 

The maxima of xT/N tot for model 2, Fig.^b), decrease with temperature as N increases. 
This indicates that there may be no finite-temperature phase transition in the thermody- 
namic limit. This conclusion is further supported in Fig. El There the location of both 
the crossing of the Binder cumulant and 0.5 and the maximum of xT/N tot is shown. Both 
estimates for the 'finite system-size' critical point agree. They both also have a reasonable 
fit, based on JV > 32, of T c = AN* with x = -0.105 and A = 2.74. 
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FIG. 5: The maximum of the susceptibility times the temperature for model 1. The slope of this 
gives an estimate for the ratio of critical exponents j/u. 

IV. POSSIBILITIES FOR SMALL WORLD MATERIALS 

Model 2, with the small-world bonds built of atomic spins, exhibits no finite temperature 
phase transition in the thermodynamic limit. Consequently, the outlook for small- world 
materials is that they will not be easy to synthesize. Note that rigorously it has been shown 
that no finite temperature phase transition can occur in one dimension. Nevertheless, there 
are a large number of effective one- dimensional materials [1] . Small- world models with fixed 
connection sizes, such as model 1, exhibit a finite temperature critical point. Building these 
small-world connections from atomic spins as in model 2 provides the small-world bonds. 
However although the small-world fixed point is unstable at finite temperatures, the flow 
from this unstable fixed point seems to be extremely slow. Consequently, the outlook for 
quasi small-world materials is better than for quasi-one-dimensional materials. 

Table I shows the effective T c e ff for a given size of material for model 2. Also shown from 
ref. jl^l are values of T C)e g for small- world connections with strengths that behave as Jir^ a , 
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FIG. 6: The summary for model 2 of the various indicators of critical behavior are shown as a 
function of the system size. Note that this fit (dashed line) suggests that a finite-temperature 
phase transition does not survive taking the thermodynamic limit. 

for different values of a. All temperatures are in units of J\. Due to the slow decrease in 
the critical point with the system size, Table I shows that even for large systems (of the size 
of meters) the system should exhibit an effective critical temperature that can still occur at 
a reasonably high temperature. 

One synthesis route for small-world nanomaterials is described here. This is a theorist's 
cartoon of the synthesis route, and consequently is not meant to be detailed. First, start 
by constructing a one- dimensional chain of atoms (or molecules), most likely embedded on 
the surface of a support material. Prepare a solution of linear molecules of all different 
lengths up to the maximum size of material to be synthesized. This could be accomplished 
by starting with linear molecules of the same length, and then having a chemical process 
that breaks the linear chains at random locations. Place a reactive agent at the ends of 
these linear segments, this agent should be reactive with the constructed one-dimensional 
chain of atoms on the surface of the material. Bring this constructed chain of atoms into 
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Material 
Size 


-^atom 


T c 

added atoms 


T 

P 

a = 0.1 


T 

f 

a = 0.2 


T 

f 

a = 0.3 


base 


1 


2.74 


1.62 


1.67 


1.69 


jim 


5 x 10 3 


1.12 


1.48 


1.31 


1.69 


mm 


5 x 10 6 


0.54 


1.37 


1.02 


0.70 


m 


5 x 10 9 


0.26 


1.26 


0.73 


0.26 



TABLE I: Extrapolated values for the finite-size effective critical temperature (for example T p , 
the location of the maximum in %T) for various sizes corresponding to the given number of atoms 
in a one-dimensional chain material. See the text for the assumptions. The last three columns 
represent the data of ref. |l2j, as explained in the text. 



contact with the solution of reactive linear molecules, forming the small-world connections. 
A further synthesis step can be performed to reduce the length of the linear molecules to 
the minimum possible distance. This could be accomplished by removing atomic segments. 
The material synthesized in this way should mimic closely the small-world connections of 
model 2. 

There are several factors that may make quasi small-world nanomaterials easier to syn- 
thesis than the synthesis route outlined above. It has been shown that not all original atoms 
(ones that are not in the small-world connections) need to be connected with small-world 
bonds for the system to be controlled by the small-world fixed point . Rather, what 

is required is just a finite density of small-world bonds. Hence, the prefactor in front of 
the power law for the number of total atoms per original atom may be made arbitrarily 
small. It is also possible to consider starting with a two-dimensional thin film, and synthe- 
sizing small-world connections on the film in a fashion similar to the one-dimensional route 
outlined above. Some studies of crossover from two- and three-dimensional fixed points to 



2lJ. However, it is anticipated that for 



small-world fixed points have been performed [20! 
any finite density of small-world connections the dominant fixed point should be the small- 
world fixed point. Furthermore, since these original one and two-dimensional systems are 
embedded in three dimensions, it is possible to make some of the longer- distance small-world 
connections very short-ranged by bending or folding the original lattice to minimize the to- 
tal length of all small-world bond connections. Finally, just as in quasi-one-dimensional 
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materials there is the possibility of allowing weak interactions in three dimensional crystals 
to stabilize the fixed point in different dimensions before the system exhibits a cross-over to 
the three-dimensional fixed point. 



V. DISCUSSION AND CONCLUSIONS 



The ferromagnetic Ising model on two small- world networks has been studied. One model, 



171 ] , has only fixed-strength in- 



related to models previously studied by others [l2 , JJ , 3, 13, 
teractions for the small-world connections. This model exhibits a finite-temperature second- 
order phase transition. We have determined from finite-size scaling of the maximum of the 
susceptibility the exponent ratio 7/1/ = 0.6(1) for the system sizes simulated. A study by 
other researchers with small-world connections that fall off as a power-law with actual dis- 
tance shows no finite-temperature phase transition, but a logarithmic decrease with Nq of 
an effective finite system size critical temperature T^eff jl^t • 

The second Ising model studied here is more realistic for the possibilities of quasi small- 
world materials. The small-world connections have been constrained to be built from spins 
using the same lattice spacing. This models what would have to be accomplished in actual 
quasi small-world materials, since atoms are the fundamental building blocks. For this 
model, no finite-temperature phase transition was found to survive the thermodynamic limit. 
In particular, in terms of the number of linear spins Nq, an effective critical temperature 
T c ,es was found to behave as Nq with x ~ —0.105. Furthermore, on general grounds the 
total number of spins scales as N tot ~ Nq for a one-dimensional system with N spins. 

From these studies a theoretical picture emerges about the difficulty in synthesizing quasi 
small-world materials. The number of atoms needed to connect small-world bonds may be 
large, and is expected to scale as a power law of the number of atoms in the none small- 
world system. However, it has been shown that not all none small-world atoms need to 
be connected with small-world bonds for the system to be controlled by the small-world 
fixed point [3, Q|- Consequently, the prefactor in front of the power law for the number of 
small-world atoms may be made arbitrarily small in principle. For strengths of small-world 
connections that decrease as a power law [12], a logarithmic fall-off of the effective critical 
temperature with system size has been found. Using some order-of-magnitude assumptions 
and the different power law results from this work and ref . , one finds that a reasonable 
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value of T Ci efr is possible for systems in the size ranges of microns, millimeters, or even meters. 
For the model studied here, where the small- world connections were built from atomic spins, 
a small power-law fall-off of the effective critical temperature was found. Thus nanoscale 
and mesoscale materials should exhibit an effective critical behavior related to the small- 
world bonds. A possible synthesis route to materials effectively governed by small-world 
fixed points was outlined. 
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